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Abstract: We analyze the dark matter problem in the context of brane cosmology. We 
investigate the impact of the non-conventional brane cosmology on the relic abundance of 
non-relativistic stable particles in high and low reheating temperature scenarios. We show 
that in case of high reheating temperature, the brane cosmology may enhance the dark 
matter relic density by many order of magnitudes and a stringent lower bound on the five 
dimensional scale is obtained. We also consider low reheating temperature scenarios with 
chemical equilibrium and non-equilibrium. We emphasize that in non-equilibrium case, the 
resulting relic density is very small. While with equilibrium, it is increased by a factor of 
0(1O 2 ) with respect to the standard thermal production. Therefore, dark matter particles 
with large cross section, which is favored by detection expirements, can be consistent with 
the recent relic density observational limits. 
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1. Introduction 

The history of the observed Universe is not known before the epoch of nucleosynthesis (i.e., 
temperature above 1 MeV). The assumption that our universe has to be homogeneous and 
isotropic is only true at late universe. Therefore, one can consider the standard cosmology 
based on Robertson- Walker (RW) metric as an effective low energy theory for a more 
fundamental one at higher energy scale. 

Warped extra dimensions have been proposed to explain the large hierarchy between 
the electroweak scale (Mw ~ 10 2 Ge V) and the Planck scale (M p i ~ 10 19 GeV). In this class 
of models, the ordinary matter is assumed to be localized on a three-dimensional subspace, 
called brane which is embedded in a larger space, called bulk. It has been emphasized that 
the brane cosmology in these models can be quite different from the standard cosmology 
of four dimensional universe. In particular, the derived Friedman equation of a brane 
embedded in five dimensional (5D) warped geometry is given by [Q] 

„ 87rG(4) / p \ k C , 

where H = a/ a is the Hubble parameter and a{t) is the scale factor, p is the energy density 
of ordinary matter on the brane while a is the brane tension. refers to the AD Newton 
coupling constant. Finally k stands for the curvature of our three spatial dimensional and 
C is a constant of integration known as dark-radiation. As can be seen from the above 
equation, H oc p rather than ^J~p as in the conventional cosmology. Thus, the evolution 
of the scale factor will be different from the standard one. This modification would be 
very relevant for the cosmological events, as dark matter relic abundance, that may occur 
during the radiation dominated phase of the early universe. 

One of the theoretical problem in modern cosmology is the existence of dark matter 
(DM). The most interesting candidate for this dark matter is a long lived or stable weakly 
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interacting massive particle (WIMP) which can remain from the earliest moments of the 
universe in sufficient number to account for the dark matter relic density. The standard 
computation of the WIMP relic density, based on the usual early universe assumptions, 
leads to the following fact: The WIMP relic density is inversely proportional with its 
annihilation cross section. However, the detection rate of this particle, which is given 
in terms of its elastic cross section with the nuclei in the detector, is proportional to its 
annihilation cross section. Therefore, in order to detect the WIMP experimenally, its cross 
section should large, of 0(1O -6 — 10 -8 )GeV -2 . Neverthless, in this case the WIMP relic 
density is quite small: Q x h 2 < 0.01 which contradicts the recent observational bounds: 
0.094 < a x h 2 < 0.128 |]. 

The aim of this paper is to study the relic density in the context of non-conventional 
brane cosmology. It is important to mention that there are some attempts for similar 
analysis in the literature [|], H, however, that work was focused on the standard thermal 
production scenario for computing the relic density. In this case, the reheating temperature 
is assumed to be much larger than the freeze-out temperature and the dark matter particles 
are in thermal equilibrium. Here, we analyze the WIMP relic density in non-conventional 
brane cosmology in case of high and low reheating temperature. We show that in the first 
scenario where the reheating temperature is high, a strong lower bound is imposed on the 
51? fundamental scale. While in case of low reheating temperature and the WIMPs reach 
a chemical equilibrium before reheating, the relic density is enhanced by one-two order of 
magnitudes. This enhancement could help to accommodate large annihilation cross section 
DM which, as mentioned, has a very small relic density in the standard scenario. As an 
example, we apply our result for the lightest super symmetric particle (LSP) which is one 
of the best candidate for cold dark matter. 

The paper is organized as follows. In section 2 we briefly review the non-conventional 
brane cosmology and point out that the consistency of the boundary conditions on the 
brane implies that the 5D metric should contain two warp factors at least, otherwise the 
relation = —p^ is identically satisfied. Section 3 is devoted for analyzing the WIMP 
relic density in brane cosmology when the reheating temperature is much higher than the 
decoupling temperature. We show that in this case the relic density is enhanced by many 
order of magnitudes which set a strong constraint on M5. In section 4 we study the DM 
relic density in brane cosmology with low reheating temperature where DM particles could 
be in chemical equilibrium or non-equilibrium. We show that the former case is the most 
interesting scenario, since the relic abundance is increased by a factor of O(10 2 ). Therefore 
it can be an interesting possiblity for accommodating the recent relic density observational 
limits with large cross sections . Finally we give our conclusions in section 5. 

2. Brane world cosmology 

In this section we discuss the effective brane cosmology from Einstein equations in five 
dimensional space-time. We consider the most general metric that preserves three dimen- 
sional rotational and translation invariance: 

ds 2 = -n 2 (t, y)dt 2 + a 2 (t, y)-f ij dx i dx j + b 2 (t, y)dy 2 , (2.1) 
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where 7^ is a maximally symmetric 3-dimensional metric with spatial curvature k = ±1,0. 
In our analysis, we identify the hypersurface defined by y = with the brane that forms 
our universe. The induced metric in this brane is the usual 4D RW metric. 

The five-dimensional Einstein equations (with bulk cosmological constant A) are given 

by 

Rab 2A ) 9ab = tf 5) f AB , (2.2) 

where Rab is the 5D Ricci tensor, R is the 5D scalar curvature and the constant Krg) is 
related to the 5D Newton's constant Gm and the 5D Plank mass Mr^\, by the relations 

4) = 8 ^ G (5) = M (5) • ( 2 - 3 ) 

Assuming an empty bulk, the energy momentum tensor is due to the matter on the 
brane, which is considered to be an infinitely thin. Thus TP is given by 

T B = M diag (_ pb; pb) pb) pb; 0) , (2.4) 

where pb and pb are the total energy density and pressure on the brane, respectively. In 
order to have a well defined geometry, the metric (i.e., the function / = n, a, b in Eq.( |2.1| )) 
must be continuous across the brane localized at y = but its derivative with respect to y 
can be discontinuous. Therefore, its second derivative with respect to y will contain Dirac 
delta function, i.e., 

/" = ?+[/']%), (2-5) 

where /" is the non-distributional part of the second derivative of / respect to y and [/'] 
is the jump in the first derivative of / across y = which is defined as 

[/']=/'(0 + )-/'(0-). (2.6) 



By matching the Dirac 5(y) function in Eq.( ^.2|) , one obtains the following relations, which 
are known as junction conditions Q: 



a Q b 
n b 



Kf 5) (2pb + 3p b ) 



(2.8) 



where the subscript stands for the evaluation at y = 0. These conditions imply that the 
5D metric should contain at least two different warp factors, otherwise we get pb = —pb 
as an identical relation which means that our universe is always vacuum dominated. This 
is not consistent and contradicts with the assumption of having matter on the brane. For 
instance, if we consider the case of embedding the AD RW metric in the following 5D 
metric: 

ds 2 = h 2 (t,y) [-dt 2 + R(t) 2 'Yi j dx i dx i ] + b 2 (t,y)dy 2 . (2.9) 
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In this case, one finds that a 2 (t,y) = R(t) 2 h 2 (t,y), n 2 (t,y) = h 2 (t,y), and the junction 
conditions (|2.7|), (|2.8[) lead to the relation 

Pb = ~Pb, (2.10) 

i.e., the vacuum energy is always dominated on the brane, which is not necessary true. It 
is worth mentioning that this result does not depend on the choice of the scale factor b(t, y) 
of the fifth dimension or on our cosmological scale factor R(i) 1 . 

Using the condition ( |2.7j ) with the components (0, 0) and (4, 4) of Einstein's equations 
in the bulk, one finds the following equation 

where C is a constant of integration. As can be seen from this equation, the Hubble 
parameter is proportional to the energy density of the brane, in contrast with the standard 
four-dimensional Friedmann equation where it depends on the square root of the energy 
density. Let us consider a brane with total energy density 

Pb = ° + P, (2.12) 

where a is a brane tension, constant in time, and p is the energy density of ordinary 
cosmological matter. This implies 

36 6 18 ^ 36^ a 2 a 4 v ' 



Then with fine tuning of brane tension, the first term in Eq.( 2.13 ) vanishes if we have 

= (2.14) 
36 6 V ' 

Furthermore, fixing the value of A in terms of and M p i as A = — QMP^/M^, one finds 

6 

which leads to the new Friedmann equation: 



87rG (4) = -^a, (2.15) 



„ 8ttG(4) / p \ k C 
tf-^U^l+i)- + (2.16) 



As can be seen from Eq.( |2.16| ), at low energies i.e., at late time, the cosmology can be 
reduced to the standard one, but in the early time the p 2 term becomes dominant, so the 
universe undergoes nonconventional cosmology. 

Before concluding this section, it is worth mentioning that the normal energy conser- 
vation equation is still valid, 

P + 3H(p + p) = 0, (2.17) 

1 For a similar observation in brane cosmology of Jordon-Brans-Dick theory, see Ref. [pi. 
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which implies that p oc a . Using this relation in Friedmann equation ( 2.16| ) one 
finds § 

ao(t) oc t q , (2.18) 

where 

* = P>>CT - (2-19) 

Therefore, within the radiation domination epoch where w = 1/3 one gets ao(t) oc in 
the non-conventional cosmology limit and ao(t) oc in the standard case limit. Thus the 
time dependence of Hubble parameter in raditation domination era is given be 

H = -, with q=\ * ■ (2.20) 

t I | /3«cr 

As mentioned in the introduction, one would expect important impact for these modifica- 
tions on the cosmological events that occurred in the early universe during the radiation 
domination era. In the next section we will analyze the effect of the non-conventional 
cosmology on the dark matter relic abundance. 

3. DM relic density in non-conventional brane cosmology 

In this section we compute the relic density of the WIMP (%) within the non-conventional 
brane cosmology. In the standard computation for the WIMP relic density, one assumes 
that x was i n thermal equilibrium with the standard model particles in the early universe 
and decoupled when it was non-relativistic. Once the x annihilation rate T x = {a x nn v)n x 
dropped below the expansion rate of the universe, T x < H, the WIMPs stop to annihilate, 
fall out of equilibrium and their relic density remains intact till now . The above {a x nn v) 
refers to thermally averaged total cross section for annihilation of XX hito lighter particles 
times the relative velocity, v. 

The relic density is then determined by the Boltzmann equation for the WIMP number 
density (n x ) and the law of entropy conservation: 

= - 3 Hn x - (of v) [{n x f - (r^) 2 ] , (3.1) 

tt = ~ 3HS > (3 ' 2) 
where n x is the WIMP equilibrium number density which, as function of temperature T, 
is given by n x = g x (m x T /2tt) 3 ^ 2 e~ mx ^ T . Here m x and g x are the mass and the number 
of degrees of freedom of the WIMP respectively. Finally, s is the entropy density. In 
the standard cosmology, the Hubble parameter H is given by H 2 = (8tt/3M 2 z ) p to be 
compared with the expression in Eq.( |2.16 ) for brane cosmology. In our analysis we will set 



k = C = in order to focus on the impact of the modification of the p dependence in H. 

Let us introduce the variable x = m x /T and define Y = n x /s with Y eq = n x /s. In 
this case, the Boltzmann equation is given by 
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In radiation domination era, the entropy, as function of the temperature, is given by s = 
^-g* s (x) rr? x x~ 3 = kix~ 3 , which is deduced from the fact that s = (p + p)/T and <?* s is 
the effective degrees of freedom for the entropy density. Therefore one finds 

ds 3s 
dx x ' 

which is the same in both cases of standard and brane cosmology. It is worth mentioning 
that since the time variation of s is given in terms of the Hubble parameter H, s(t), 
unlike s(T), is expected to be different in brane cosmology from the usual expression in 
the standard case. This indicates that the time-temperature relation is modified in brane 
cosmology. As known in the standard cosmology this relation is given by 



1 \^ M ll 9 

4 « = 2\/^ r <3 ' 5) 

In brane cosmology with dominant p 2 term in the Hubble parameter H, the time-temperature 
relation takes the following form 




=M T-*. (3.6) 

This means that within the brane cosmology, the cooling of the universe becomes much 
slower than it in standard cosmology. In the standard case, the following expression for 
the Boltzmann equation for the WIMP number density is obtained 

f-jfcK-^-ny, (3.7) 



where H(x) is given by H = y 4 \®^z x x 2 = y/k^x 2 and is the effective degrees of 

freedom for the energy density. Therefore, one obtains the following expression for the 
Boltzmann equation for the WIMP number density (g* s ~ g* is assumed): 

In brane cosmology the Hubble parameter is given by H = (fox -4 + A^x -8 ) 1 / 2 where 
&3 = 7r 4 (jr 2 m 8 < /(32400M|). Thus, the Boltzmann equation in brane cosmology takes the 
form 

(f ), - -^§ m " <»* 4 + tr /2 •> < y2 - • < 3 ' 9 » 

It is worth noticing that in the limit of k$ — » (i.e., a — * oo), the above equation 
tends to the standard Boltzmann equation in Eq. (|3.8j ). Therefore, at early times, the 
universe undergoes a nonstandard brane cosmology till it reaches a temperature, known 
as transition temperature Tt where the universe sustains the standard cosmology. This 
transition temperature is defined as Q 

p(T t ) = 2a => T t = 0.51 x 10~ 9 m| GeV (3.10) 
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In order to analyze the brane cosmology effect on the WIMP relic density, one should 
assume that the freeze out temperature of the WIMP (Tp) is higher than the transition 
temperature, i.e., Tp > Tt. Therefore, one finds 



M 5 < 1.57 x 10° I — 

Xp 



771., 



2/3 



(3.11) 



Since the WIMPs freeze out at temperature Tp <C m x , they are non-relativistic and there- 
fore the averaged annihilation cross section can be expanded as follows: 



(oT n v) 



66 

a H , 

x 



(3.12) 



where a describes the s-wave annihilation and 6 comes from both s- and p- wave annihi- 
lation. To obtain the present WIMP abundance Y^, we should integrate the Boltzmann 
equation for the WIMP number density from xp (the decoupling temperature) to Xoo — oo 
(present temperature). It is important to notice that this integral must be divided to two 
parts from xf to Xt where the non-conventional brane cosmology is applied and from xt to 
oo where the universe undergoes the standard cosmology. Thus, one obtains 



Y. 



-i 

ocb 



71 '9* m 



x F 



( a + — ) I r- + z 

x \k 2 



-1/2 



dx + 



•ft 



a 66 

X A x^ 



(3.13) 



Here we have used the usual assumption that Y eq <C Y and Y XF 3> Y^. Evaluating the 
above integrals, one finds 



Y. 



oob 



45 



M p i m x 



a ( - 2 Fi 

x 



3t/^6|sinh-' 



-k 3 



115 

4' 2' 4' k 2 x i 




sinh 




(3.14) 



where 2F1 [a, 6, c, z] is the Hypergeometric function, which is a solution of the hypergeomet- 
ric differential equation: z(l — z)y" + [c — (a + 6 + l)z]y — aby = 0. As can be seen from 
Eq. ( p. 14 ) that for xt = xf the expression of Y^ 6 coincides with the standard known result 



for Y oo\i namely ^ 
given by 



n x h 2 



Px 



The relic abundance of the WIMP is 



2.9 x 10 8 K, 



771^ 



(3.15) 



-GeV, 

where the critical density p c is given by p c ~ 10 _5 /i 2 GeVcm -3 and h is the Hubble 
constant, h ~ 0.7. Furthermore, the p x is defined as p x = tii x so^oo where so — 2900 cm -3 
is the present entropy density. As in the standard scenario, the relic density of the WIMP 
is inversely proportional to its annihilation cross section. However, unlike the standard 
case, it depends explicitly on WIMP mass since k 3 /k 2 oc m x . 

The freeze-out temperature can be determined from the freeze-out condition: 



A(x F ) = cY?{x F ), 



(3.16) 
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where A is given by A = Y x — Y x q and c is a constant of order unity. In standard cosmology, 
this condition implies that 

0.0765 c m x M p ig x (a + 6b/x F ) , , 

xp = in , (3.17) 

\/xf9*{xf) 

which can be solved iteratively to determine the value of xf- It turns out that for m x ~ 
0(100) GeV, xf ~ 0(25). In brane cosmology, one can easily show that xf is obtained by 
iterative solution of 

xf _ ln 0-0765 c m x M p ig x x 3 / 2 (a + Qb/x F ) ^ ^ 

In this case, one finds that xp is smaller than the above value obtained within the standard 
cosmology. Also, it turns out that xf is sensitive to the scale M5. For example with 
M 5 ~ 10 6 one gets x F ~ 0(7). 

Let us introduce the factor R = (£l x h 2 ) / (£l x h 2 ) s that measures the enhancement/suppression 
in the relic abundance due to the brane cosmology. From Eq.( 3.15| ), one finds 

a _i_ 36 

R= XF ^ 



1 a 1 So 
-r — -t- -7 



'£b ( sinh" 1 ( A /£^ 2 ) - sinh- 1 (Jt^U+aih.F, 

(3.19) 

This ratio could be larger or smaller than one depending on the values of the annihilation 
cross section parameters a and b and also on the masses m x and M5. In order to ana- 
lyze the impact of the non-conventional brane cosmology on the relic density result, we 
consider, as an example, the lightest supersymmetric particle (LSP), which is one of the 
most interesting candidates for the WIMP. As is well known, in most of the parameter 
space of the supersymmetric models the LSP is mainly pure Bino. Therefore, it is mainly 
annihilated into lepton pairs through t— channel exchange of right-handed sleptons. The 
p-wave dominant cross section is given by @ 



6~ 87ra 2_L 1 (3.20) 
m x v 1 + x i R ) 

where xj r = m~ /m x and a\ is the coupling constant for the U{l)y interaction. Thus, for 
m x ~ m l R ~ 100 GeV, one finds b ~ 0(1O~ 8 ) GeV -2 , which in the standard cosmology 
scenario leads to Q. x h 2 > 0.1. 

As advocated above, in brane cosmology the relic density {Sl x h 2 ) is quite sensitive 
to the value of the fundamental scale M5 which should satisfy the upper bound given in 
Eq.fl3~Tl|). Therefore, with m x ~ 0(100) GeV and x F ^ 0(10), one finds 



M 5 < 10 7 . (3.21) 

Furthermore, the fact that the transition process from non-conventional cosmology to con- 
vention cosmology should take place above the nucleosynthesis era (i.e., Tt > 1 MeV) 
impose the following lower bound on M5: 

M 5 > 1.2 x 10 4 . (3.22) 
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Log[M 5 ] 

Figure 1: The enhancement /suppression factor R = (fl x h 2 )b/ '(fl x h 2 ) s as a function of the five 
dimensional scale M 5 (GeV) for m x = 100 (solid curve), 200 (dashed curve) and 500 GeV (dotted 
curve). 

In Fig. |l], we present the prediction for the factor R as a function of the scale of the 
five dimensions, M5, for different values of m x , namely we consider m x = 100, 200 and 500 
GeV. As can be seen from this figure, for M5 < 10 6 the brane cosmology effect is quite large 
and the factor R becomes much larger than one. In this case the resulting relic density 
(Q x h 2 )b may exceed the WMAP results (at 95% confidence level) j| 

n x h 2 = 0.1126±g;gi|i. (3.23) 

Moreover for M5 > 5 x 10 6 , the ratio R becomes less than one and a small suppression 
for (£2 x h 2 ) s can be obtained. This brane enhancement or suppression for the dark matter 
relic density could be favored or disfavored based on the value of the relic abundance in the 
standard scenario. If (Q x h 2 ) s is already larger than the observational limit, as in the case of 
bino-like particle, then a suppression effect would be favored and hence M5 is constrained 
to be larger than 5 x 10 6 GeV. However, for wino- or Higgsino-like particle where the 
standard computation usually leads to very small relic density, the enhancement effect will 
be favored and the constraint on M5 can be relaxed a bit ||. In general, it is remarkable 
that in this scenario the dark matter relic density imposes a stringent constraint on the 
fundamental scale M5. 

4. DM relic density in brane cosmology with low reheating temperature 

In the standard computation for the DM relic density that we have adopted in the previous 
section, it was assumed that the reheating temperature Tjuj is much higher than the WIMP 
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freeze-out temperature i.e., T R h 3> Tp. In this case, the reheating epoch has no impact 
on the final result of the relic density. However, it is well known that the only constraint 
on Trh is Trh ^ 1 MeV in order not to spoil the successful predictions of the big bang 
neucleosynthesis. Therefore, in principle it is possible to have a cosmological scenario with 
low reheating such that Trh < Tp. In this case the predictions of the relic abundance of 
the WIMP are modified as emphasized in Ref.|7]] and recently in Ref.S. 

As in the previous section and to emphasize the effect of the brane cosmology, we 
assume that the WIMP freeze-out temperature is larger than the transition temperature 
which is also larger than the reheating temperature, i.e., Tp > Tt > T R h- Within a low 
reheating temperature scenario, the relic density depends on whether WIMPs are never in 
chemical equilibrium either before or after reheating or they reach chemical equilibrium 
but they freeze-out before the completion of the reheat process. The resultant relic density 
in these two scenarios are quite different and are also different from the one derived by the 
standard computation with a large reheating temperature. 

Let us start our analysis with the case of non-equilibrium production and freeze-out. In 
this case, at early times the number density n x is much smaller than n x and the Boltzmann 
equation (3.7) can be written as 



0.02095 ±± —{al nn v) x 6 e~ Ax . (4.1) 
dx V 9*s J Hx x 

Although this equation is valid only at early times, it can be approximately integrated in 
the full range of x, namely from x = to x = oo due to the exponential suppression in 
the right hand side. Thus for the standard cosmology (where H oc x~ 2 ), one can easily 
integrate this equation and obtains 

Yoo S = 0.02095^ ~^glg- 3/2 M pl m x (| + 3&) . (4.2) 

The Yqo is related to the mass density of x particle today as follows, at reheating we have 
P x ( t Rh) = m x n x (T RH ) = 2*. g^ s {TjiH)rn x Y 00 T^ H . After the reheating the universe is 
radiation dominated and the following relation is satisfied M: 

P X ( T now) _ T RH p x (T RH ) ^ ^ 

Pr(Trh)' 

Therefore, in this case Q x h 2 is proportional to the annihilation cross section instead of 
being inversely proportional as in case of high reheating temperature. 

Now we consider this scenario of low reheating with non-equilibrium production and 
freeze-out in brane cosmology. The Boltzmann equation is still given by Eq.([4.1|), but with 
H = (k 2 x~ 4 + /C3X- 8 ) 1 / 2 in the range of x G [0, xt] and with the usual Hubble parameter 
H = y/hzx -2 between xt and x = oo. Integrating this equation one finds 

Ycob ~ 0.02095 x 10~ 6 ^J^9l9* 3/2 M pi m x ( 9 - 46 a + 37 - 8 b ) ■ ( 4 - 4 ) 

Here we have used m x ~ 100 GeV and M5 ~ 10 6 GeV, as an example, to do the integration 
numerically. However, we have checked the result of for different values of m x and M5. 
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It turns out that is diminished significantly for M5 < 10 6 . As can be observed from 
equations ( |4,2| ) and (4^4) that this non-equilibrium scenario produces a very suppressed 



relic density, particularly in brane cosmology. Furthermore, the assumption that n x <C n x q 
impose a sever constraint on the annihilation cross section. Thus, one can conclude that 
within this scenario, it is not possible to account for the dark matter experimental results. 

Now we turn to the second scenario in which the annihilation cross section of the 
WIMP is large and hence it reaches the chemical equilibrium before reheating. In this 
case, the computation of the relic density Q x h 2 is very close to the standard case with high 
reheating temperature. At the early times i.e., when T > Tp, the WIMP's are very close 
to equilibrium. So, as usual, one can use the variable A(x) = Y(x) — Y eq {x) to write the 
Boltzmann equation as 



A' 



f(x)A (2Y eq + A) 



(4.5) 



where f(x) is given by f(x) 
obtains 



(Y eq ) 



661 



A ~ 



. Then by neglecting A' and A 2 , one 



(4.6) 



2f(x)Y eq 

At late time, when T < Tp one gets Y 3> Y eq and hence we can use the approximation 
Y 2 — (Y eq ) 2 ~ Y 2 in Eq,([3.8|) and integrate from Tp down to Trh to determine Y(Trh). 
In the standard cosmology, one finds 



1 



Y(x RH ) 



1 



Y(xp) 



Kg* 
45 



M p im x 



+ 



3b 



X p 



(4.7) 



If it is assumed that there is no entropy production for T < Trh, then there is no WIMP 
production for temperature below the reheating temperature. Thus, the present value of 
Y is given by Y(xrh) up to an overall correction due to the fact that the reheating process 
is not complete at Trh 0- Here, two comments are in order: i) As mentioned above, 
Y(xf)^ 1 ~ Y eq (xp)^ 1 which is of order O(10 9 ), so its contribution to Y{xrh)~ 1 in Eq.(^7|) 
can be neglected respected to the second term, ii) Since Trh < Tf (i.e., xrh > xf), one 
can approximate Eq^^) and finds that the relic abundance VL x h 2 is very close to the one 



obtained by using the standard calculation with high reheating temperature, namely 

-1 



n^h 2 ~ 1.1 x la 



36 

, + 

\Xf x\ 



n I a 6b 



(4. 



which implies that unless annihilation cross sections are quite small ( < 10~ 8 ), one gets, 
as usual, very small relic density. 

In brane cosmology, Eq.(|3.S|) describes the correct Boltzmann equation that should be 
used. Also, as in the previous scenario, one has to integrate this equation from xp to xt 
using brane cosmology feature, and from x t to xrh using the standard cosmology feature. 
In this respect, one finds 



1 



1 



Y(xrh) b Y(x F ) 



71 9 * Ail 





sinh \ \ — 



sinh 



3^-2 
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One can easily check that the three parts of the second term in the above equation are of 
the same order and give the dominant contribution to Y(xrh)~ 1 ■ in this case, one can 
show that the relic density (for m x = 100 GeV and M5 = 10 6 GeV) is given by 

a x h 2 ~ 1.1 x 10~ 7 (95.2 a - 4.12 by 1 . (4.10) 

From this equation it can be easily seen that even with large annihilation cross section 
0(1O~ 6 — 10~ 8 ), we are able to obtain the cosmologically interesting value S7 x /i 2 ~ 0.1. 
This implies that the scenario of brane cosmology with low reheating termperature and 
chemical equilibrium WIMPs is the most interesting model for dark matter. It provide an 
interesting possibility for having dark matter with large cross section (hence their detection 
would be possible in future DM experiments) with suitable relic aboundanc. 

5. Conclusions 

In this paper we have analyzed the relic abundance of cold dark matter in brane cosmology. 
We have pointed out that the consistency of the boundary conditions that account for the 
presence of the brane implies that the 5D metric should contain, at least, two warp factors. 
In case of one warp factor, we have emphasized that the relation = —p^ is identically 
satisfied, which contradicts the possibility of having matter on the brane. 

We have also studied the impact of brane cosmology on the cold dark matter relic 
density. We investigated this effect in two different scenarios, namely when the reheating 
temperature is higher or lower than the freeze-out temperature. We showed that with 
high reheating temperature, the relic density is enhanced with many order of magnitude 
for M5 < 10 6 . This imposes one of the strongest constraints on the scale of large extra 
dimensions. In case of low reheating temperature, we have considered the possibility that 
WIMPs are in chemical equilibrium or non-equilibrium, which depends on the value of their 
annihilation cross section. We emphasized that if WIMPs are in chemical non-equilibrium, 
then their relic density is very small and they can not account for the observational limits. 
While in case WIMPs reach chemical equilibrium before reheating, we showed that the relic 
density is enhanced by two order of magnitudes than the standard thermal scenario result. 
This enhancement can be considered as an interesting possibility for accommodating dark 
matter with large cross section, which is favored by the detection rate experiments. 
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